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ABSTRACT 

We discuss the fragmentation of a filamentary cloud on the basis of a 
1-dimensional hydrodynamical simulation of a self-gravitating gas cloud. The 
simulation shows that dense cores are produced with a semi-regular interval in 
space and time from one edge to the other. At the initial stage the gas near 
one of the edges is attracted inwards by gravity and the accumulation of the 
gas makes a dense core near the edge. When the dense core grows in mass up 
to a certain amount, it gathers gas from the other direction. Accordingly the 
dense core becomes isolated from the main cloud and the parent filamentary 
cloud has a new edge. This cycle repeats and the fragmentation process 
propagates towards the other edge. The propagation speed is a few tens of 
percent larger than the sound speed. According to the theory, the age difference 
for the northwest-most and southeast-most cores in TMC-1 is estimated to be 
0.68 pc/0.6 km s _1 = 10 6 y. The estimated age difference is consistent with that 
obtained from the chemical chronology. 

Subject headings: hydrodyanamics, instabilities, ISM: individual objects: Taurus 
Molecular Cloud 1. 
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1. INTRODUCTION 



The star formation process seems to propagate by several mechanisms on both large 
and small scales (see, e.g., a review by Elmegreen 1992). A classical example of propagating 
star formation is the linear sequences of stellar subgroups in Orion OB associations (piaauw 



|1964|) . A larger scale example is the super giant shell in the Large Magellanic Cloud 



( |Wcsterlund fc Mathewson 1966| ). The shell consists of HI clouds and bright blue stars, and 
surrounds two supernova remnants and stellar associations. Similar shells are also seen in 
our Galaxy. In these examples, massive stars trigger star formation in the neighborhood 
through supernova explosion and UV radiation. 

The star formation process may also propagate in a dark cloud. Recently Hirahara 
et al. (1992) have found that dense cores are successively formed. According to them, a 
filamentary cloud, Taurus Molecular Cloud 1 (TMC-1) is elongated in the NW-SE direction 
and contains five dense cores called A-E from NW to SE. The cores have different chemical 
ages and are older in NW than in SE. The age difference is of the order of 10 6 y. This 
indicates propagation of dense core formation. 

If cores A-E are formed successively, it is likely that a causal mechanism works there. 
However it could be neither supernova explosion nor UV radiation from young stars; no 
influential stars are found in the vicinity of TMC-1. Although various models are proposed 
for sequential star formation on large and small scales, most of them invoke young massive 
stars. They cannot be applied to TMC-1. 

In this paper we propose a new mechanism for sequential core (star) formation 
which operates in a filamentary molecular cloud. In our model, formation of a dense core 
triggers that of another core. Once a core is formed, the gravitational field is changed and 
accordingly the gas flows. Another core is produced by the gas flow and the cycle repeats. 



-4- 



Using a 1-dimensional (ID) hydrodynamical model we demonstrate how this mechanism 
works. 

In section 2 we construct a ID model for the dynamics of a filamentary cloud. Some 
characteristics of the ID model are investigated analytically. It is shown that the ID 
model can explain the fragmentation of a filamentary cloud due to the self-gravity of the 
gas. Numerical simulations of the ID model are shown in section 3. They demonstrate 
that cores are produced sequentially with semi-regular temporal and spatial intervals in a 
filamentary cloud. The apparent propagation speed of fragmentation is somewhat larger 
than the effective sound speed. In section 4 we summarize observations so far reported for 
TMC-1 and apply our model to TMC-1. The apparent propagation speed of fragmentation 
in TMC-1 is of the order of the turbulent velocity. Since the filament length and velocity 
dispersion are 0.68 pc and 0.6 km s -1 in TMC-1, the time scale for the propagation of 
fragmentation is estimated to be 10 6 yr. This estimate is consistent with the age of cores 
calculated from their chemical abundances. In section 5 we discuss the dynamical and 
chemical evolution of a filamentary cloud during fragmentation. 



When a filamentary cloud fragments, the gas motion is mainly in the direction parallel 
to the filament axis. Thus, we approximate the gas motion to be 1-dimensional. The 
equation of motion and the equation of continuity are expressed as 



2. 



1-DIMENSIONAL MODEL 



2.1. Model Equations 




(1) 
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and 



respectively, where v, p, P, and ib are the velocity along the filament axis, the density, the 
pressure, and the gravitational potential, respectively. The Poisson equation is arranged to 
be 

d 2 ib ib , , 

where G and i? are the gravitational constant and a length scale for the filament radius. 
The derivative in the radial direction in the Poisson equation is replaced with a constant 
coefficient [ the second term in the left hand side of equation (3) ]. Equations (1) thorough 
(3) denote the gas motion in a filamentary cloud under the 1-dimensional approximation. 



2.2. The Filamentary Cloud Model and Its Stability 

In this subsection we discuss a model for a infinitely long (d/dz = 0) filamentary cloud 
and its stability. A steady state solution of equations (1), (2), and (3) is expressed as 

p(z) = po = const. (4) 

v(z) = , (5) 

and 

ib(z) = -AttGpqR 2 . (6) 

We consider a small perturbation around this equilibrium. The density perturbation is 
assumed to be in the form of 



p(z, t) = po + Pi exp (ikz + iut) 



(7) 
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Then the dispersion relation is expressed as 

- 2 = -4vrGp or ^, + c s V. (8) 

Equation (8) is qualitatively similar to the approximate dispersion relation for a filamentary 
cloud flHanawa et al. 1993| ). 

^ = -^G Pc -^— + c 2 k 2 , (9) 
1 + kti 

where the equilibrium density distribution is assumed to be 

po = p c (l + r 2 /8H 2 )- 2 . (10) 
This similarity justifies our ID approximation and suggests R ~ H. 
In our ID approximation, the model is unstable whenever 

AnGp R 2 > c 2 . (11) 

We think that inequality (11) is satisfied in a gravitationally bound filamentary cloud. For a 
2-dimensional model of a filamentary cloud there is a condition of magnetohydrodynamical 
equilibrium in the radial direction; 

B 2 

4ttG Pc H 2 = c 2 + 2n 2 H 2 + — — (1 + cos 2 9) , (12) 

167rp c 

where Q c , B c and 9 denote the angular velocity at the filament center, the magnetic field 
density at the center, and the pitch angle of the magnetic fields at r = oo (see Hanawa 
et al. 1993 for further details). Equation (12) means that the gravitational force should 
be equated with the sum of the thermal and magnetic pressures and the centrifugal force. 
If a filamentary cloud is permeated with magnetic fields and/or rotates around the axis, 
inequality (11) is fulfilled. For example, TMC-1 has velocity gradient along the minor axis 
which can be ascribed to rotation around the major axis (e.g., Dlano, Walmslcy fc Wilsori 



|1988| ). In the following we assume inequality (11) except when otherwise stated. 
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The growth rate and wave number of the most unstable mode are derived to be 

(-iu) m ax = y/foGpo - C 8 /R, (13) 
and 



k - * 
' R \ 



RV^Upi _ j (14) 



C s 

from equation (8). The values of (— iu) max , k max , and A = 2ir/k max are summarized in 
Table 1. Equations (13) and (14) are used to interpret the ID nonlinear simulation of 
fragmentation. 



2.3. Periodically Distributed Cores 



Equations (1) - (3) have a more general steady state solution where the density changes 
periodically. In hydrostatic equilibrium the density is expressed as 

1>(z) - #))■ 



p(z) = p(0)exp{- 



-}• 



(15) 



Substituting equation (15) into equation (3) we obtain 



= 47rGp(0)exp{- 



dz 2 tt" ' y c s 

In the following we assume that the density has the minimum value, p(0) = p min at 
z — 0, and hence the gravitational potential is the highest at z — 0, with ^(0) = ip n 
Integrating equation (16) multiplied by dip/dz we obtain 



(16) 



dip 2 



|^(*)|2 ^(0) 



-87rG p(0) c s 2 {exp(- 



ij(z) - ^(0) 



:i7) 
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Since dip/dz is expressed as a function of ip, we can obtain z as a function of ip for given 
p(0) and ^(0) which satisfy 

- 47rGp(0) J R 2 < ^(0) < . (18) 
When ip(0) = — 4irGp(0)R 2 , the density and gravitational potentials are spatially constant. 



NUMERICAL SIMULATION 



3.1. Numerical Methods 



The differential equations (1) - (3) are integrated numerically in Lagrangian coordinates 

with a fully implicit scheme. We introduce a Lagrangian coordinate, s, defined by 

ds = p(z) dz . (19) 

Using s we rewrite equations (1) - (3) into 

(I). = «• (*» 

O, - -(£)-•• ™ 

P = c s 2 p, (22) 

P - *(£)", (23) 



>ds 

dg p c il) 



ds ' pR 2 



+ 4ttG Pc , (24) 



j- s = {Pol 9)9- (25) 

We followed the time evolution of a Lagrangian mesh point, Zi (0 < % < n) using the 
difference equations, 

Zj = zf A) + btv j , (26) 
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vj = ^ (old) + A*(-c 2 Pj+1/2 7 J '- 1/2 - ft-) , (27) 

Pc t-^S 

where 

= Pc(^^) _1 (28) 

^ +1/2 - — ^ As 

- 47rGp J+1/2J R 2 , (29) 

3 2p c As 1 ; 

We take the fixed boundary condition, 

Zq = const, and z n = const. , (31) 

at the both ends. The number of the Lagrangian mesh points is n = 1.1 x 10 4 in most of 
numerical simulations shown in this paper. The time step is taken to be At = 5.0 x 10 -3 . 

In the following we take units of R = 1, p = 1, and 2-nGpo = 1 except when 
otherwise stated. 



3.2. Standard Model (Model 1) 

As a typical example we describe model 1 where the sound speed is taken to be 
c s = 1.0. The initial gas distribution is assumed to be 

p = 0.1 p (for - 50 < ,2 < 0) (32) 
= 1.0 po (for < z < 50) (33) 
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and 



v = . 



(34) 



The time evolution of model 1 is shown in figure 1 with grey scale representation. The 
abscissa and ordinate are the z-coordinate and time, respectively. 

At an early stage with t ~ 2, the gas near the left edge of the cloud is attracted 
by the cloud gravity and flows to right. A dense core is produced near z ~ 4 at t ~ 4. 
The maximum density at t = 6 is p = 3.574 at z — 5.433. The core is in a quasi-static 
equilibrium after t > 8. As the core grows in mass, it attracts gas also from its left side. 
The first core becomes isolated from the main cloud at t ~ 10 and the cloud has a new 
edge. 

The second core is produced in a similar way at t ~ 12. Cores are produced 
sequentially with a time interval of At ~ 8 and spatial interval of Az ~ 9. The spatial 
interval is nearly equal to the wavelength of the fastest growing mode in this model 
(^max = 9.763). The apparent propagation speed of fragmentation, Vf rag = Az/At ~ 1.1, 
is a little faster than the sound speed. Nakamura, Hanawa & Nakano (1993) estimated the 
propagation speed to be 



from a linear stability analysis (see their Appendix 4). Substituting equations (13) and (14) 
into equation (35) we obtain 



Equation (36) gives a good estimate, Vf rag = 1.23 c s for model 1. 

Figure 2 shows the density and velocity distribution at t — 20. The density profile 
has peaks at z — 9.00, 15.73, and 23.30, where the peak density is 18.06, 9.19, and 2.60, 
respectively. The first and second cores are in a quasi-static state and the velocity is almost 




(35) 




(36) 
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constant in each core. These cores accrete gas from both sides and standing shock waves 
are formed at the core surfaces. The formation of the third core is in progress. The first and 
second cores attracts each other and merge into a single core at a later stage. This means 
that the periodic distribution of gas discussed in subsection 2.3 is also unstable against 
merging. 

The exact epoch of merging depends on the details of the initial density distribution 
and even on the number of the mesh points. Although a trial computation with fewer mesh 
points is quite similar to model 1 in the early stage, the first and second cores merge earlier 
in the trial computation. This is because fewer mesh points give a larger amplitude of 
numerical noise. Compare model 1 with models 4 and 5 for the dependence on the initial 
density distribution. 



3.3. Dependence on c s 

In order to study the dependence on c s we constructed models 2 and 3 where the 
sound speed is taken to be c s = 1.1 and 1.2, respectively. The initial density distribution 
of models 2 and 3 is the same as that of model 1. 

The time evolution of models 2 and 3 are shown in figures 2 and 3, respectively. Models 
2 and 3 are qualitatively similar to model 1. The spatial interval between cores is Az ~ 10 
and ~ 12 for models 2 and 3, respectively. The propagation speed of fragmentation is 
Vf r ag c± 1.2 and 1.4 for models 2 and 3, respectively. When c s is larger, the cores are wider 
in separation and fragmentation propagates faster. These relations qualitatively agree with 
equations (14) and (36). 
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3.4. Dependence on Initial Density Distribution 

We constructed models 4 through 7 to study the dependence on the initial density 
distribution. The initial density profile has a sharp edge at z = in models 1, 2, and 3. 
The sharp edge is introduced for making a model which is as simple as possible and unlikely 
to be realized in interstellar space. In models 4 through 7 the initial density increases from 
0.1 to 1.0 with a constant density gradient (dp/dz) in the region —a < z < a where a is 
taken to be a = 2.5, 5, 20, and 50 in models 4, 5, 6, and 7, respectively. The sound speed is 
taken to be c s = 1.0 in models 4 through 7. Figures 5, 6, 7, and 8 show the time evolution 
of models 4, 5, 6, and 7, respectively. Sequential fragmentation also takes place in models 4 
through 7. 

Model 4 is not much different from model 1. A major difference between them is the 
distance between the first and second cores at t = 24. It is shorter in model 1 than in 
model 4. A small difference in an early stage is amplified by merging instability. 

Model 5 is also similar to models 1 and 4 except that formation of the first core is 
delayed by At ~ 2 and shifted to the right by Az ~ 3 in model 5. The delay is due to the 
low density gradient at the initial stage and the shift is because the effective edge is shifted 
to the right in model 5. Model 5 is qualitatively similar also in the formation of the second 
and third cores to models 1 and 4. A sharp edge at the initial stage is not essential for the 
propagation of fragmentation. 

Model 6 shows also the propagation of fragmentaion. In model 6 cores are formed both 
to the left and right of the first core. Propagation to the left is slower than that to the 
right. This is because the mean density is lower and the self-gravity is weaker in the left of 
of the first core. Propagation to the right terminates at t ~ 28 when it reaches to z = 50, 
the boundary of computation. Propagation to the left will terminate after it reaches to a 
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less dense region (see also figure 6). 

In model 7 the initial density increases from 0.1 to I in the region —50 < z < 50. 
The first core is formed at z = 50 where the initial is the maximum. Since the reflection 
boundary is assumed at z — 50, framentation propagates both to the left and right from 
z = 50 although only propagation to the left is shown in figure 8. 

As seen models 4 through 7 propagation of fragmentation is ubiquitous in the 
gravitational instability of a long filamentary cloud. Once a core is formed, subsequet cores 
are formed in general on the left and right sides. When the first core is formed near the 
edge of a cloud, propagation to the diffuse side diminishes (see figure 6). The first core may 
often be formed near an edge of a filamentary cloud since the large density gradient thereof 
triggers formation of the first core. In such a case cores are formed from an edge to the 
other in a long filamentary cloud. 



3.5. Comparison with 2D Simulation 

Bastien (1983) followed the non-linear fragmentation of a non-rotating non- magnetized 
cylindrical cloud with a 2D hydrodynamical code for the first time. The numerical 
simulation was extended by Rouleau & Bastien (1990), Bastien et al. (1992), and Arcoragi 
et al. (1992). They followed the collapse and fragmentation of an initially uniform 
cylindrical cloud having density, po, diameter, D, and length L. 

The qualitative features of their simulation depend mainly on J , the ratio of the 
gravitational to the thermal energies. When J > J c , the cloud is not massive enough 
to be gravitationally bound. When J > J sp i n di C ) the cloud collapses as a whole into a 
thin needle and does not fragment. The case of J c < Jo < J sp mdie corresponds to our 
ID simulation which implicitly assumes hydrostatic balance in the radial direction. When 
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J c < J < Jspindio two condensations are formed on the axis. They are symmetric with 
respect to the cloud center and move toward the center while their density increases. They 
collide with each other at the center when J is close to J c . The whole cloud becomes a 
thin needle when J is relatively large. 

Our simulation is qualitatively similar to their simulations up to the formation of 
the two condensations, one of which is the mirror image of the other. A condensation is 
produced near the cloud edge in both simulations as far as we know from the published 
figures. Further evolution is, however, very different between their simulations and ours. 
Second generation condensations are not formed in their simulation. The difference comes 
mainly from the fact that the hydrodynamical equilibrium in the radial direction can be 
realized only with difficulty for an isothermal non-magnetized gas. Only when the mass per 
unit length has a certain value (= 2c s 2 /G) can an isothermal non-magnetized cylindrical 
cloud be in equilibrium. The difference may in part be due to the numerical scheme. It is 
difficult for an explicit code to follow further evolution after the formation of high density 
cores. Low spatial resolution may induce the coalescence of high density condensations at 
an earlier stage. We think that condensations are formed successively on the axis in a 2D 
numerical simulation if the initial cloud is in a quasi-static equilibrium and the numerical 
code can follow the long term evolution. 



4. APPLICATION TO TMC-1 

In this section we briefly summarize the observations of TMC-1 with emphasis on its 
star formation history. 

Taurus Molecular Cloud 1 (TMC-1) is one of filamentary clouds in the Taurus dark 
cloud complex, whose distance from the Sun is about 140 pc (Elias 1978). The filamentary 
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structure of this cloud was first recognized by Little and his collaborators (1978) through 
mapping observations of HC5N. Since then TMC-1 has been extensively studied using 
various molecular lines such as NH3, CS, HC3N and HC7N (Tolle et al. 1981; Snell, Langer, 
& Frerking 1982; |01ano, Walmsley fc Wilson 1988j ). It is now established that the filament 



is extended from SE to NW with an apparent size of 17' x 2', which corresponds to a 
linear size of 0.68 x 0.08 pc ( [Hirahara et al. 1992| ; |(Jlano, Walmsley fc Wilson 1988| ). A 



significant velocity gradient (3.8 km s _1 pc -1 ) is observed along the minor axis (Olano et 
al. 1988). Such a velocity structure is interpreted as possible rotation or overlapping of 
several filaments. The global magnetic field in this region was studied by Moneti et al. 
(1984). They carried out optical and infrared polarimetry in this region, and found that the 
magnetic field direction is from SW to NE, being perpendicular to the filament. 

This source has also been a good target for detailed studies on chemistry in dark 
clouds, since a number of molecular species are detected. Little et al. (1979) first pointed 
out that a large chemical gradient is seen along the major axis of the filament. They 
observed the NH3 and HC3N lines, and found that NH3 is intense in the NW part whereas 
HC3N is intense in the SE part. However the origin of such a chemical gradient has long 
been a puzzle. 

Hirahara et al. (1992) mapped TMC-1 using the CCS lines with a spatial resolution 
of 40". Since the optical depths of the spectral lines of CCS are not very high, these lines 
are suitable for detailed studies on the structure and chemistry of the cloud. They reported 
that TMC-1 is not a uniform filament but a chain of several dense cores spaced almost 
regularly (the left panel of Figure 9). Five cores are identified, each of which has a diameter 
of 0.04 - 0.10 pc. These cores can also be seen in the NH3 map (the right panel of Figure 
9). They measured the cloud density accurately by observing the optically thin emission of 
the C 34 S ( J = 1 — and J = 2 — 1) lines, and found that the density tends to increase from 
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SE to NW. The density averaged over the 40" beam is 4 x 10 4 cm" 3 at the cyanopolyyne 
peak (Core D), whereas it is 4 x 10 5 cm -3 at the NH 3 peak (Core B). The size and density 
of each core are comparable to typical values for dense cores reported by Benson & Myers 
(1989). 

From these results, the chemical gradient seen in TMC-1 is ascribed to a systematic 
chemical difference among the dense cores. Hirahara et al. (1992) suggested that the 
chemical difference originates from the core age. According to gas-phase chemical model 
calculations ( Millar fc Hcrbst 1990| ; Suzuki et al. 1992), carbon-chain molecules like CCS 



are abundant in the early stages of cloud evolution, while NH3 is only abundant in the 
late stages. On the basis of this model, the NW cores are considered to be older than the 
SE cores. In fact, the density tends to be higher in the NW core. Furthermore, the IRAS 
source (04381+2540) exists near the NW core, indicating that star-formation has already 
taken place. 

The chemical lifetime of carbon-chain molecules is about 10 6 yr according to chemical 
model calculations (Leung, Herbst & Huebner 1984; Tarafdar et al. 1985; Suzuki et al. 
1992). This time scale is mainly related to the conversion time scale from C to CO. 
When a cloud is in a diffuse stage with A v < 2, the major forms of carbon are C + and 
C, because molecules are mostly photodissociated by the interstellar ultraviolet radiation. 
As the cloud becomes opaque due to gravitational contraction, photodissociation becomes 
ineffective, and hence, C + and C are consumed to form a stable molecule, CO, through 
various chemical reactions. The conversion time scale from C to CO in the opaque region 
is inversely proportional to the cosmic ray ionization rate and almost independent of the 
cloud density. Because a high abundance of C + and C is essential to the production of 
carbon-chain molecules, carbon-chain molecules can only exist in the early stages of cloud 
evolution. The lifetime of carbon-chain molecules is almost comparable to the production 
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time scale of NH 3 . Therefore, the abundance ratio between carbon-chain molecules and 
NH 3 is a good indicator of evolutionary stages during the 10 6 yr after the cloud becomes 
opaque. 

On the basis of the above scenario on chemical evolution, the age difference between 
the cyanopolyyne peak position and the NH 3 peak position in TMC-1 is expected to be 
about 10 6 yr, i.e., the cores in TMC-1 are successively formed within this time scale. On the 
other hand, our model simulation suggests that the propagation speed of core formation is 
slightly larger than the effective sound speed. The velocity width observed in TMC-1 is 0.4 
- 0.8 km s _1 , depending on the position. When we adopt 0.6 km s^ 1 as the effective sound 
speed, the difference in ages between the cyanopolyyne peak and NH 3 peak is estimated 
to be 7 x 10 5 yr. This value is almost consistent with the time scale estimated from the 
chemical differences. These results further support the prediction that the cores in TMC-1 
are successively formed from NW to SE. 



5. PHYSICAL AND CHEMICAL EVOLUTION OF CORES 

Thus far, TMC-1 is unique as a filamentary cloud indicating sequential core formation. 
Our mechanism for the sequential core formation, however, seems to work in a filamentary 
molecular cloud. The uniqueness of TMC-1 may be solely due to the fact that TMC-1 has 
been observed extensively in various molecular emission lines. In this section we discuss 
guidelines for searching for filamentary clouds which may contain cores with different ages. 

As discussed in the previous section, the chemical evolution begins as a result of the UV 
shielding due to fragmentation and contraction. Let us suppose that the visual extinction 
is A v ~ 3 at the onset of the chemical evolution. This means that the penetration of UV 
radiation to the core is reduced by a factor of a thousand at the onset. Since the visual 
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extinction is well correlated with the integrated intensity of a molecular emission line, we 
can estimate the brightness of a core at the onset of the chemical evolution. Cernicharo & 
Guelin (1987) evaluated the correlation between the visual extinction and the integrated 
intensities of 12 CO, 13 CO, and C 18 to be 

J( 12 CO) = 5.0 ±0.5- (A v -0.5±0.2)Kkms" 1 , (37) 

J( 13 CO) = 1.4 ±0.2- (A v -0.7±0.3)Kkms" 1 , (38) 

and 

/(C 18 0) = 0.28 ± 0.5 • (A v - 1.5 ± 0.3) Klaus" 1 , (39) 

for the Taurus region. The core of A v = 3 is already opaque for 12 CO, semi-transparent for 
13 CO and transparent for C 18 0. This means that the decay of CCS and production of NH 3 
take place in a dense core seen in the C 18 map. 

A core which can be identified in a C 18 map is thought to have the a density ^ 3x 10 4 
cm -3 . The dynamical time scale of such a dense core is short, 

1 

Tdyn ~ V^P 

= I.55xl0 5 f ^ ;Y 1/2 yr- (40) 

V3xl0 4 cm- 3 ^ 3 v ; 

Core collapse and formation of a subsequent core takes ~ 5 and ~ 10 dynamical time scales, 
respectively. Since the time scale for chemical evolution is ~ 10 6 yr, it is comparable to 
time scales for collapse and propagation of fragmentation only when riu 2 ^ 3 x 10 4 cm~ 3 . 
Accordingly we can observe cores of different chemical compositions only in a dense core 
which can be selected from a C 18 map. 

In fact, such a chemical difference is observed in several clouds in the Taurus region. 
Good examples are TMC-1C and TMC-1A. These clouds are visible on the C 18 map 
(Cernicharo & Guelin 1987). TMC-1C has an elongated structure from SE to NW with 
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a linear size of 0.3 x 0.1 pc flCcrnicharo, Guclin fc Asknc 1984| Yamamoto et al. 1993Q . 



The NH 3 map reported by Benson & Myers (1989) has a peak at the southeast part of the 
core, whereas the CCS map has an additional peak at the northwest part ([yamamoto et al" 



1993|) . Thus a significant chemical gradient can be seen along the major axis of the cloud. 
In this case, the SE core is considered to be older than the NW core, and hence, the cores 
are being formed successively from SE to NW. 

A similar chemical gradient is also seen in the TMC-1A region, as reported by 
Cernicharo et al. (1984). The distribution of NH 3 has a peak at the infrared source position 
(IRAS 04365+2535). On the other hand, the emission from carbon chain molecules, HC3N 
and HC5N, tends to be intense in the 5' north position, although these molecules are not 
fully mapped in this region. Therefore, the southern part of TMC-1A seems to be in an 
advanced stage of cloud evolution. It is interesting to note that a high-velocity outflow 
has recently been detected toward the infrared source ([Terebey, Vogel fc Myers 1989| ; 



Hirano 1993|) , indicating that star formation has already taken place in the southern region. 



Although a few examples of a chemical gradient in the dense cores have been reported, as 
mentioned above, more extensive observations of NH 3 and carbon-chain molecules would 
be indispensable for observational confirmation of the sequential core (star) formation in 
filamentary clouds proposed in the present paper. 
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c s 




A 




v frag 


0.9 


0.7559 


8.312 


0.5142 


0.756 


c s 


1.0 


0.6436 


9.763 


0.4142 


1.231 


c s 


1.1 


0.5345 


11.756 


0.3142 


1.512 


c s 


1.2 


0.4225 


14.871 


0.2142 


1.711 


C s 


1.4 


0.1008 


62.358 


0.0142 


1.985 


C s 


1.4142 
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Table 1: Dispersion Relation for ID Model Filamentary Cloud 
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Fig. 1. — The time evolution of model 1. The abscissa and ordinate are the z-coordinate 
and time, respectively. The density is denoted by the darkness. The grey scale is shown on 
the top of the panel. The sound speed is taken to be c s = 1.0. 

Fig. 2. — Density and velocity distribution in model 1 at t — 20. The upper and lower 
panels show the density on a logarithmic scale and velocity in a linear scale, respectively. 
The abscissa is the z-coordinate. 

Fig. 3. — The same as figure 1 but for model 2 where the sound speed is taken to be c s = 1.1. 

Fig. 4. — The same as figure 1 but for model 3 where the sound speed is taken to be c s = 1.2. 

Fig. 5. — The same as figure 1 but for model 4 where the initial density density increases 
from 0.1 to 1.0 in the region —2.5 < z < 2.5. The sound speed is taken to be c s = 1.0. 

Fig. 6. — The same as figure 1 but for model 5 where the initial density density increases 
from 0.1 to 1.0 in the region — 5 < z < 5. The sound speed is taken to be c s = 1.0. 

Fig. 7. — The same as figure 1 but for model 6 where the initial density density increases 
from 0.1 to 1.0 in the region —20 < z < 20. The sound speed is taken to be c s = 1.0. 

Fig. 8. — The same as figure 1 but for model 4 where the initial density density increases 
from 0.1 to 1.0 in the region —50 < z < 50. The sound speed is taken to be c s = 1.0. 

Fig. 9. — Total integrated intensity maps of CCS (left) and NH 3 (right) emission toward 
TMC-1 (Hirahara et al. 1992). The lowest contour level and the contour interval are 
0.4 Kkms" 1 . Five cores (A - E) are identified in the filament. 



